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ABSTRACT 



Context. V4332 Sgr experienced an outburst in 1994 whose observational characteristics in many respects resemble those of the 
eruption of V838 Mon in 2002. It has been proposed that these objects erupted because of a stellar-merger event. 
Aims. Our aim is to derive, from observational data, information on the present (10-15 yrs after the outburst) nature and structure of 
the object. 

Methods. We present and analyse a high-resolution (R ss 21 000) spectrum of V4332 Sgr obtained with the Subaru Telescope in 
June 2009. Various components (stellar-like continuum, atomic emission lines, molecular bands in emission) in the spectrum are 
analysed and discussed. We also investigate a global spectral energy distribution (SED) of the object mostly derived from broadband 
optical and infrared photometry. 

Results. The observed continuum resembles that of an ~M6 giant. The emission features (atomic and molecular) are most probably 
produced by radiative pumping. The observed strengths of the emission features strongly suggest that we only observe a small part 
of the radiation of the main object responsible for pumping the emission features. An infrared component seen in the observed SED, 
which can be roughly approximated by two blackbodies of ~950 and ~200 K, is ~50 times brighter than the M6 stellar component 
seen in the optical. This further supports the idea that the main object is mostly obscured for us. 

Conclusions. The main object in V4332 Sgr, an ~M6 (super)giant, is surrounded by a circumstellar disc, which is seen almost edge- 
on so the central star is obscured. The observed M6 spectrum probably results from scattering the central star spectrum on dust grains 
at the outer edge of the disc. 
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1. Introduction 

V4332 Sgr is an un usual object, whose nova-like eruption was 
observed in 1994 (Martini et alj 1 19991) . Discovered as a K- 
type (super)giant, it showed a rapid cooling on a time scale 
of months and declined as a very cool M-type (super)giant. 
This evolution suggests that V4332 Sgr has a similar na- 
ture to V838 Mon, wh ich erupted in 2002. As discussed by 
Tylenda & Soker (2006), neither the classical nova mechanism 
nor the He-shell flash model can explain the observed outbursts 
of these objects. At present the stellar collision-merger scenario 
proposed in Soker & Tylenda (2003) and further developed in 
Tylenda & Soker ( 2006) is the most promising hypothesis for 
explaining the nature of these eruptions. 

Although it is now more than a decade after the outburst 
of V4332 Sgr, it still shows very unusual observational char- 
acteristics. In the optical and near-IR, the object displays a 
very unique emission-line spectrum superimposed on a weak 
continuum resembling a photospheric spectrum of a cool star 
(Baneriee et al. 2003; Baneriee & Ashok 2004; Tvlendaetal. 
2005; Ki meswengeril2006l) . The emission spectrum shows very 
low excitation and is composed of lines from neutral elements 
(e.g. KI, Nal, Cal) and molecular bands of TiO, ScO, AlO. 
Excitation temperatures derived from analyses of the emission 
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* Based in part on data collected at Subaru Telescope, which is oper- 
ated by the National Astronomical Observatory of Japan. 



spectrum range in 200-800 K. The origin of this spectrum is not 
clear. 

V4332 Sgr is very bright in the infrared. Spectroscopic stud- 
ies showed strong absorption features of water ice at 3.05 and 
6 /im, of silicates and alumina at ~10 /mi, and possibly of 
methane ice at 7.7 //m, as well as the fundamental band of CO 
in emission at 4.67 /jm (Baneriee et al. 2004, 2007). Early stud- 
ies of the spectral energy distribution showed that, apart from 
the stellar-like component dominating the observed continuum 
in the optical, the object displays a bright infrared component 
with a characteristic temperature of ~750 K and a luminosity at 
least 15 tim es high er than that of the stellar one dTvlenda et al.l 
2005; Ba neriee & A shok 2004). A possible circumstellar disc or 
ring orbiting the central object was discussed in several stud- 
ies (Baneriee & Ashok 2004; Baneriee et al. 2004; Tvlenda et al. 
2005). V4332 Sgr remained relatively constant in the optical 
brightness between 2002 and 2005. In 2006-2007, the object 
slowly faded by ~1.5 mag in the B, V, and R bands, remaining 
almost unchanged in the / bandj. 

In this paper we further investigate the nature of V4332 Sgr. 
Our study is based on new high-dispersion spectrosopic obser- 
vations in June 2009 with the Subaru Telescope. We also discuss 
the global spectral energy distribution (SED) of the object de- 
rived from different measurements done in the optical, near-IR, 
and far-IR. We show that the emission spectrum is most likely 
produced by radiative pumping. From our analysis of the opti- 
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cal photospheric spectrum, emission features, and the SED, we 
conclude that the central object in V4332 Sgr, most probably a 
cool (super)giant, is hidden in a circumstellar disc seen almost 
edge-on. 



2. Observations 

2.1. Optical spectroscopy with Subaru/HDS 

High-resolution spectroscopy of V4332 Sgr was obtained 
on 20 09 June 16 with th e High Dispersion Spectrograph 
(HDS, iNoguchi et all |2002|) on the 8.2-m Subaru Telescope. 
Observations were performed in the standard StdRb setting with 
two different echelle angles, resulting in a combined spectrum 
covering the range 5287-8070 A at a nominal resolution of 
R =* 21 000. Owing to the inter-chip gap and some irregulari- 
ties of the instrument's CCDs, the spectrum is missing regions 
6640-6708 A and 7340-7360 A, and has a few other minor gaps. 
Five individual exposures were obtained with a total time of 100 
min. 

Data were reduced with IRAF0 using standard procedures 
for echelle spectroscopy (IWillmarth & Bar nes 1994). All the 
spectra were wavelength-calibrated using a Th-Ar spectral lamp 
and flux-calibrated using our obse rvations of HR 7596 and spec- 
trophotometric fluxes from Hamuy et al. (1994, 1992). We as- 
sess an overall accuracy of the relative flux calibration at ~ 15% 
and that of the wavelength calibration at 0.01 A. 



2.2. NIR photometry 

Photometric JHKs observations of V4332 Sgr were carried 
out on 2009 May 18 with the CAIN-3 camera installed on 
the 1.52 m Carlos Sanchez Telescope (Observatorio del Teide). 
Observations were performed in a sequence of exposures 4x10 
s, 4x7 s, and 12x3 s for the JHKs filters, respectively. Data 
were reduced with the CAIN data reduction script^] developed 
in IRAF by the Instituto de Astrofisica de Canarias. The data 
were calibrated to standard magnitudes by using photometry 
from the All-Sky Point Source Catalog of the Two Micron All 
Sky Survey (2MASS). The averaged magnitudes for V4332 Sgr 
are /=12.60±0.16, #=10.79+0.15, and tf s =9.12+0.15. 



2.3. Submillimeter observations of CO(3-2) 

V4332 Sgr was also observed in the 12 CO(3-2) (345.796 
GHz) line on 2009 May 14 using the James Clark Maxwell 
Telescopqj (under t he programme S09AI01). The HARP array 
dBuckle et al.ll2009h was used in the conventional stare mode. 
The telescope's half-power beamwidth at the observed fre- 
quency is about 15". The observing method was beam-switching 
with a chop throw of 60" in azimuth. The ACSIS autocorrela- 
tor was used as a back-end providing a spectral resolution of 
0.42 km s . The total on-source integration time was 25 min 



2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 

3 http://www.iac.es/telescopes/cain/reduc/caindr.html 

4 The James Clerk Maxwell Telescope is operated by The Joint 
Astronomy Centre on behalf of the Science and Technology Facilities 
Council of the United Kingdom, the Netherlands Organisation for 
Scientific Research, and the National Research Council of Canada. 



and the typical system temperature during observations was of 
T sys = 402 K. 

We did not detect any emission in the LSR velocity range 
between -500 km s" 1 and 400 km s" 1 for all of the HARP re- 
ceptors. The upper limit on the CO(3-2) emission at the position 
of V4332 Sgr is 3cr = 78.7 mK (per 0.42 km s" 1 channel) on 
the antenna's temperature scale. This corresponds to 141 mK on 
the scale of the main beam temperature (T m b), if the main beam 
efficiency of T] m \, = 0.56 is assumed. 



3. The optical spectrum of V4332 Sgr in 2009 

The Subaru/HDS spectrum of V4332 Sgr appears similar to 
spectra that have b een reported for the post-outburst evolution of 
this object to date dBanerjee & Ashold2004lTvlenda et all2005b 



iKime s wen geif 2006 ) . Strong lines of neutral atoms and molecu 
lar bands, all seen in emission, are displayed on a cold M-type 
photospheric spectrum. The best part of our spectrum is shown 
in Fig.Q] Because the object is very faint, the raw spectrum has 
a very poor signal-to-noise ratio (S/N) in the original resolution; 
the average S/N per pixel in the continuum is of about 0.7 in 
the 5500 A region and ~7 near 8000 A. To identify photospheric 
features, the spectrum was smoothed applying boxcar smoothing 
with a box-size of 21 pixels. Fortunately, most of the emission 
features are very strong, and they can be studied in detail even 
at the original resolution. The photospheric spectrum, and the 
atomic and molecular emission features, as well as interstellar 
features, are all presented in detail in the following sections. 

3. 1 . Photospheric spectrum 

A photospheric spectrum, characteristic of a late type star, can be 
easily recognized in Fig . Q] When compared with earlier spec- 
troscopic observations dBaneriee & Ashokll2004_ iTvlenda et alJ 
2005; Kimeswenger 2006), the present spectrum displays much 
more pronounced wide and deep absorption bands in the red 
part. This obviously suggests a later spectral type and a lower 
effective temperature. We constrained the spectral type of the 
present photospheric spectrum. The Subaru/HDS spectrum, cor- 
rected for interstellar extinction with Eg-v - 0.32 (see Sect.|4]i, 
was compared to a grid of observed spect ra of late type gian ts 
wit h luminosity classes III— II, taken from Bagnulo et al. (2003 Jj 
and lGunn & Strvkendl983h . The result is shown in Fig. [2 where 
the giant (luminosity class III) spectra of types M4, M5, M6, and 
M7 are plotted against the spectrum of V4332 Sgr. From this 
comparison we can classify the V4332 Sgr spectrum as M5-6. 
This can be compared with s pectral types K8-M3 derived from 
observational data in 2003 dTvlenda et al.M 2005: Kimeswenger] 
2006). The decline in the spectral type most probably took place 
in 2006-2007 when the object was gradually fading in BVR, as 
noted in Sect. Q] The derived spectral type of the 2009 spec- 
trum corresponds to an effective temperature of T e g =* 3300 K 
(Richichi et aljfl999l) . We have also fitted standard photomet- 
ric spectra of luminosity class III to the BVR C I C J photometric 
data presented in Sect. [7] (We used the same procedure as in 
Tvle nda et al.| [2005) The best fit gives a spectral type of M6.2 
and r eff =* 3200 K. 

We were not able to constrain the luminosity class of the 
object, except that it is a giant. There are no clear signs of ab- 
sorption features of CaH, characteristic of spectra of late-type 
dwarfs. 



see also http://www.sc.eso.org/santiago/uvespop/field_stars_uptonow.html 
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Fig. 1. The best part of the spectrum of V4332 Sgr obtained with Subaru/HDS in June 2009 (red). The spectrum has been smoothed 
and dereddened with Eb-v = 0.32 and Ry = 3.1 (see text for details). The strongest emission features are identified. Parts affected 
by telluric absorption bands are indicated by horizontal bars at the bottom. 




Fig. 2. The 2009 spectrum of V4332 Sgr (red - the same as 
in Fig. ED) compare d to spe ctra of M-type gian ts taken from 
iBagnuloet ail d2003l) 5 and lGunn & Strvkerl d 1983b . The regions 
suspected of being affected by some "extra" molecular absorp- 
tion are indicated by vertical green marks. 



The observed spectrum of V4332 Sgr differs from any stan- 
dard spectrum of a late-type giant. In addition to the obvi- 



ous presence of emission features, we noticed extra absorption 
bands not observed (or being much weaker) in stars of spectral 
types close to M6 III. Most striking is the B-X absorption band 
of VO in the range ~7334-7534 A, which is more character- 
istic of giants earlier than M6. Some extra absorption can also 
be found in the range ~6985-7050 A; it can be assigned to the 
satellite branch of TiO y (0,0) (cf. lKaminskietaL 2009), but this 
identification is very tentative. The regions contaminated by the 
extra absorption features are indicated in Fig. [2] 

A possible explanation for the difference between the ob- 
served spectrum of V4332 Sgr and the standard spectra of gi- 
ants could be the peculiarities in the object's metallicity. We 
compared the Subaru/HD S spectrum with a grid of MARCS 
spectra dGustafsson et al . 2008) with T e $ - 3200 K, logg = 
0.0, M = 1 Mq, standard chemical composition, and metallic- 
ity in the range -2 < [Fe/H] < +0.5 (see the table on 
http://marcs.astro.uu.se). The shape of the observed spectrum 
is best represented by the model spectrum with [Fe/H] = -1.5. 
The fit is still not perfect but no extra absorption is needed in 
the 7500 A region; however. iMartini et alJ (1999) obtained fits 
to the spectra observed during the 1994 outburst using models 
with solar abundances. Our suggestion of a low metallicity in 
V4332 Sgr is very tentative. The object is not a typical giant, so 
peculiarities in the spectrum are expected. Besides, as discussed 
in Sect. [7] we probably do not observe the object directly, but 
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rather a part of its spectrum that is scattered on a disc matter. 
Scattering can affect the shape of the spectrum. 

We did not find any atomic absorption features that could be 
identified as photospheric lines. A photospheric absorption line 
of Ca 1 6572 A, relatively strong in late type giants, may be the 
cause of the irregular shape of the observed Ca I emission profile, 
as discussed in Sect. 13.21 No photospheric signatures were found 
within the profiles of the Na I and K I doublets. 

The lack of atomic absorption lines in the photospheric spec- 
trum does not allow us to derive a reliable value of the radial 
velocity of the object. As can be seen from Fig. [1] most of 
the photospheric absorption bands coincide with narrow com- 
ponents in emission. However, in a few cases relatively clean 
band-heads without emission component can be seen and their 
position measured. This is the case for the following TiO y 
bands: (3,2) R x (6849+1.0 A), (4,3) Ri (6917.6+0.5 A), (5,4) 
Ri (6987.4+ J. A), (0,0) s R 2 i (7036.7+0.5 A), and (3,4) R 3 
(7818.6+0. 5 A). These measurements give rough estimates of 
the (heliocentric) radial velocity of the object, i.e., -40, -56, 
-45, -67, and -60 km s , respectively. These values imply 
a heliocentric radial velocity of the source of the photospheric 
spectrum of -56 ±16 km s . 

3.2. Atomic emission lines 

The strongest atomic emission lines seen in the Subaru/HDS 
spectrum are listed in TableQ] They are characterized in terms 
of observed wavelength, heliocentric radial velocity (V/,), full 
width at half-maximum (FWHM), integrated flux, and equiva- 
lent width (EW). For the KI and Nal doublets, we also present 
'cleaned' fluxes, which are fluxes measured for the lines after re- 
moving all atmospheric and interstellar absorption features from 
their profiles. The fluxes are not corrected for interstellar red- 
dening. We do not present any measurements of EW for the 
Na I doublet since the observed continuum in the region of these 
lines is very faint. Some of the strongest emission lines are in- 
dicated in the spectrogram presented in Fig.Q] Their profiles are 
displayed in more detail in Fig. [3] 

The profiles of the emission lines are triangular in shape and 
their cores can be reproduced well by a Lorentzian profile (see 
Fig.[3j middle panel); however, the observed wings are signifi- 
cantly less extended than in the Lorentzian profile. The emission 
lines appear to be fairly symmetric. Only with the mirror profile- 
inverting procedure or by over-plotting a theoretical profile on 
the observed ones can some asymmetry be noticed (see Fig. [3]). 
The Ca I line shows a clear sign of asymmetry, with a stronger 
blueshifted half of the profile, but this may be an effect of a non- 
flatness of the underlying continuum or a contamination by a 
photospheric absorption line (if present). 

3.3. Molecular bands in emission 

All the identified molecular emission features are listed in 
Table |2 Columns (l)-(2) of the table contain the observed and 
laboratory wavelengths of the feature. Columns (3)-(4) show the 
name of the molecule and identification of the electronic system, 
vibrational band, and the branch forming the head. Column (5) 
presents the heliocentric radial velocities of those features for 
which a unique assignment was reliable. Column (6) shows the 
integrated fluxes of the emission lines (for which reliable mea- 
surements were possible), while column (7) lists their equivalent 
widths. Comments on the nature of the features are given in col- 




V h (km s-) 



Fig. 3. Profiles of the strong resonance lines in the spectrum of 
V4332 Sgr. Spectra are rescaled. The vertical lines mark the cen- 
tre of the K I lines at -69 km s -1 and the centre of interstellar ab- 
sorption features of Nal at -7.2 km s _1 . Top: The Nal profiles, 
clipped for clarity. The top spectrum shows the intensity-inverted 
profile of the HI 21 cm line (see text). Middle: The observed K I 
lines and a theoretical Lorentzian profile. The strongest telluric 
absorption lines are indicated. Bottom: Ca I and Rb I lines, spec- 
tra slightly smoothed to reduce the noise. 



umn (8), while the sources of the laboratory wavelengths can be 
found in column (9). 

The laboratory wavelengths refer to the band-head positions. 
When a sharp band-head is recognized in the spectrum, the ob- 
served wavelength corresponds to the position of the band-head. 
These values have been used when determining the radial veloc- 
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Table 1. Measurements for the resonance emission lines seen in the Subaru spectrum of V4332 Sgr. 



-Ihb (air) 


'lobs 


Ion 


v h 


FWHM" 


FWHM" 


Flux 


Flux clean 6 


EW 


[A] 


[A] 




[kms- 1 ] 


[A] 


[kms- 1 ] 


[erg s -1 ever 2 ] 


[erg s~ l cm" 2 ] 


[A] 


5889.950 


5888.82+0.14 


NalD, 


-57.5±7.2 


2.72+0.05 


138.5 


(1.81±0.01)e-15^ 


(2.21±0.06)e-15 c 




5895.924 


5894.77±0.07 


NalDi 


-58.7±3.6 


2.68+0.06 


136.3 


(2.00±0.01)e-15 c 


(2.22±0.07)e-15 c 




6572.779 


6571.12±0.02 


Cal 


-75.7+0.7 


1.21+0.13 


55.2 


(7.34±0.43)e-16 




9.5+1.3 


7664.911 


7663.13±0.03 


KI 


-69.7+1.1 


3.54+0.02 


138.5 


(5.73±0.02)e-14 


(6.00±0.02)e-14<' 


285.4+33.1 


7698.974 


7697.22±0.01 


KI 


-68.3+2.9 


3.36+0.04 


130.9 


(4.89±0.01)e-14 rf 


(4.83±0.02)e-14 d 


267.2+17.7 


7800.268 


7798.52±0.01 


Rbl 


-67.2+0.2 


1.58+0.31 


60.7 


(1.49±0.08)e-15 




6.0+0.6 


7947.603 


7945.98+0.03 


Rbl 


-61.2+1.3 


1.01±0.08 rf 


38.1 s 


(3.54+0. 13)e-16" 




0.85+0.1" 



" from a fit of a Lorentzian profile; full width at half-maximum is not corrected for the instrumental profile (FWHMs 
* integrated flux with interstellar and strongest telluric features removed from the profile 



0.15 A) 



results from a deblending procedure with a double Lorentzian profile; simple integral for both doublet lines gives summary flux of 



(3.26±0.06)e-15 erg s" 1 cnT 2 and cleaned flux of (3.48±0.06)e-15 erg s" 1 cm 
d line blends with a TiO band; width and flux measurements uncertain 



-2 



ities. In the other cases, the observed wavelengths are the po- 
sitions of apparent maxima of the emissions. In the case of the 
VO B 4 II-X 4 X~ bands, the four emission features correspond to 
the four sub-bands formed by the main branches concentrated 
around the sub-band origin. 

The emission fluxes were measured by integrating the sur- 
face above the continuum. The continuum position has usually 
been determined by an inspection of the VLT/UVES spectra of 
the M5 III, M6 III, and M7 I II stars (HP 214952 , HD 189124, 
and HD 18242, respectively) (Bag nulo et al.ll2.Q03b . In rare cases 
of gaps in the comparison spectra, a theoretical spectrum was 
synthesized to exclude the possibility of strong underlying ab- 
sorptions. Large errors in the flux measurements of the VO B 4 LT- 
X 4 2T (0,0) are mainly due to uncertainties in the underlying con- 
tinuum. 

A detailed analysis of the profiles of the bands is presented 
in Sect. |5] 

3.4. Interstellar features 

The only features that can be ascribed to the interstellar medium 
(ISM) are the Na I 5889/5895 and KI 7664/7698 doublets in ab- 
sorption. No diffuse interstellar bands were found in the spec- 
trum. We determined radial velocities, FWHMs, and EWs of the 
interstellar features by fitting Gaussian profiles. Before perform- 
ing a fitting procedure, the appropriate regions of the spectrum 
were normalized with high order polynomials to correct for the 
broad emission features that form a steep baseline for the inter- 
stellar features. The results of Gaussian fits are listed in Table[3] 
The errors of the EWs reflect only a statistical error, which was 
determined from a series of independent measurements for each 
line. The FWHMs given in Table [3] are not corrected for the in- 
strumental broadening. The FWHM of the instrumental profile, 
as measured from the narrowest telluric lines seen in the spec- 
trum, is on average 8.0 km s , so the profiles of the interstellar 
features are mostly resolved. 

As can be seen in Table [3] there is a shift of about 4 km s _1 
between the peaks of the Nal and KI lines. The shift may be 
physical, but it can also be partially caused by the fact that the 
Na I absorption lines are saturated or close to saturation. 

The interslellar absorption features seen in the spectrum of 
V4332 Sgr can be compared to a profile of the H I 21 cm emis- 
sion line in the direc tion of the star. From the 21 cm LAB Survey 
(Kal berla et al.l 120051) . we extracted a spectrum for a position 
nearest to the position of V4332 Sgr, i.e. at / = 13!50, b = -9?50 



Table 3. Results of Gaussian fits to the ISM absorption features 
found in the Subaru/HDS spectrum of V4332 Sgr from June 
2009. 



^obs 


Ion 


v b 


FWHM 


EW 


[A] 




[km s _1 ] 


[kms- 1 ] 


[A] 


5889.950 


NaID 2 


-7.13 


15.78 


0.40 (±0.04) 


5895.924 


NalD, 


-7.32 


13.73 


0.34 (±0.02) 


7664.911 


KI 


-3.95" 


11.34° 


0.18 (±0.03)" 


7698.974 


KI 


-3.27 


9.73 


0.12 (±0.01) 



feature blends with a telluric line; measurements uncertain 



(V4332 Sgr: I = 13!63, b = -9!40). The half-power beam-width 
of the LAB Survey is 36', so the position of V4332 Sgr is lo- 
cated close to the centre of the beam for the extracted position. 
The profile of the H I line is composed of several narrow compo- 
nents, tightly blending into one complex profile. The intensity- 
inverted and rescaled spectrum of HI is shown in Fig. [3] (top 
panel). A Gaussian fit to the overall profile gives a central ve- 
locity of V/, = -6 km s . The profiles of the absorption lines 
of K I and Na I lie well within the profile of H I, The interstel- 
lar absorption lines seen in the optical spectrum indicate that the 
line of sight towards V4332 Sgr crosses most, if not all, of the 
ISM in this direction. This conclusion agrees wit h an analysis of 
extinction of V4332 Sgr in Kimes wenged (120071) . 



4. Reddening, distance, radial velocity 

Analyses of spectroscpic and photometric observations of 
V4332 Sgr during its outburst in 1994 have le d Martini et al.l 
dl999l) . lTvlenda et all (l2005h . and iKimeswengen (120061 ) to con- 
clude that the reddenning to the object is Eb-v = 0.32 + 0.02, 
E B -v = 0.32 ± 0.10, and E B . V = 0.37 ± 0.07, respectively. 

The reddenning can be estimated from a relation between 
EWs of atomic inte rstellar absorption lines a nd Eb-v- Using the 
relation derived in Munari & Zwitter (1997), our measurement 
for the K I 7699 line from Table [3]results in E B -v = 0.45 + 0.05. 
We did not use the Na I D lines since they reach zero flux in our 
spectrum and are probably saturated. 

We also used the Galactic Dust Extinction Service which 
estimates Eb-v from the 100 fim dust emission mapped by 
IRAS and COBE/DIRBE. This tool is based on a method that 



http://irsa.ipac.caltech.edu/applications/DUST/ 
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Table 2. Molecular emissions 



[A] 



•!|ab 

[A] 



Molecule Identification 



Velocity Flux 

[kms- 1 ] [10~ 16 erg s~' cm" 2 ] 



EW Comment 

[A] 



Reference 



5449 


5448.233 


TiO 


a(0,l)R 2 








2.8+0.5 


38+22 


noisy 




5474 


5469.3 


VO 


C 4 2r-X 4 E 


- (1,0) R 4 




0.6±0.3 


5+4 


weak 




5736 


5736.703 


VO 


c 4 x--x 4 s 


" (0,0) R 4 




< 1.0 




weak 




5795 


5794.6° 


CrO 


B 5 n-x 5 n(i,o)Ri 




-70+20* 


7.0+0.6 


28+5 




2 


6036 


6036.17 


ScO 


A 2 n 3/2 -X 2 E + (0,0) 


«QlG 


-75±10 b 


3.9±0.6 


14.5+4.5 






6058 


6051.8'' 


CrO 


B 5 n~x 5 n (0,0) r. 






3.5+1.0 


12.5+3.5 


noisy 


2 


6078 


6079.30 


ScO 


A 2 n 1/2 -X 2 E + (0,0) 


^0 


-75±10 


7.0+0.6 


21.0+5.0 






6133 


6132.097 


YO 


A 2 ni /2 -X 2 E + (0,0) 


k 2 , 




1.2+0.6 


2.5+0.7 






6161 


6158.52 


TiO 


r'(P.0)*i 








4.5+0.8 


31+8 






6190 


6186.32 


TiO 


Y (o,o) r 2 








1.1+0.8 


13+7 


noisy 




6212 


6214.93 


TiO 


Y (0,0) R 3 








0.8+0.8 


5+5 






6352 


6351.29 


TiO 


7 (2,0) R, 








<0.7 








6400 


6394.2 rf 


CrO 


B 5 n-x 5 n (o,i) Ri 




-70±20 b 


4.8+0.7 


10+6 




2 


6418 


6414.742 


TiO 


y(3,i)R, 












uncertain 




6449 


6451.7" 


CrO 


B 5 n-x 5 n (i,2) r, 






<1.0 




uncertain 


5 


6478 


6478.953 


TiO 


7 (4,2) R, 








<1.5 




weak 




6717 


6714.477 


TiO 


y(i,0)Ri 








2.7+0.5 


14+4 






6772 


6772.3'' 


CrO 


B 5 n-x 5 n (0,2) r. 










blend with TiO 




6785 


6781.815 


TiO 


r(2,DRi 








3.1+0.5 


8+3 


blend with CrO 




6837 


6836.5« 


CrO 


B 5 n-x 5 n(i,3)Ri 










uncertain 




7052.5 


7054.256 


TiO 


r (0,0) r 3 






-80 b 


6.6+1.5 


14+4 






7086.2 


7087.566 


TiO 


7 (0,0) R 2 






-80 b 


5.3+0.9 


14+5 






7129.4 


7125.510 


TiO 


7 (0,0) R, 






-80 b 


17.6+1.2 


44+4 


blend (1,1) R 3 




7165 


7158.850 


TiO 


y(U)Ra 








1.3+0.4 


3+1 


noisy 




7199.3 


7196.373 


TiO 


7 (2,2) R 3 








2.8+0.7 


5+1 


blend (1,1) R, 


3 


7269.3 


7269.984 


TiO 


7 (2,2) Ri 








<1.1 




weak 


1 


7381 




VO 


B 4 n-x 4 s- 


(1,0) F 2 - 


-F 2 




10.0+0.7 


5+2 






7413 




VO 


B 4 n-x 4 E- 


(1,0) F 3 - 


-F 3 




10.5+1.0 


4+1.5 






7590.4 


7589.285 


TiO 


r(o,i)R 3 








6.5+3.0 


6.5+3 


telluric Oi 


3 


7740.7 


7743.045 


TiO 


7 (2,3) R 3 






-91+39 






noisy 


3 


7826.3 


7828.153 


TiO 


y(2,3)Ri 






-71±11 






noisy 


4 


7864 




VO 


B 4 n~x 4 s- 


(0,0) Fi- 


-F, 


-75±5 b 


26+2 


12+1 






7907 




VO 


B 4 n-X 4 E" 


CCO) F 2 - 


-F, 


-75±5 b 


27+8 


11 + 1 






7946 




VO 


B 4 n-x 4 i~ 


(0,0) F 3 - 


-F, 


-75±5 b 


36+8 


15+4 


blend with Rb I 




7987 




VO 


B 4 n-x 4 s- 


(0,0) F 4 - 


-F 4 


-75±5 b 


29+10 


11+4 







References: (1) see references in lKaminski et all d2009h : (2) lRoseiJd 19701) : (3) ISchwenkeldl998l) : (4) lPhillipsld"l973l) : (5) lHocking et alUl98Ch 



" remaining four heads are at positions 5796.0 (R 2 ), 5797.5 (R 3 ), 5801.0 R,, and 5806.8 R 5 dRosenll 19701) 

b values determined from simulations of the rotational contour 

c remaining four heads are at positions 6053.3 (R 2 ), 6054.8 (R 3 ), 6058.5 (Rt), 6063.5 (R 5 ) ( Rosedll970h 

d remaining four heads are at positions 6396.2 (R 2 ), 6397.8 (R 3 ), 6401.4 (R4), 6407.7 (R 5 ) dRoseril 19701) 

e wavelengths of the band heads derived from spectroscopic constants of Hocki ng et all ( fl980h : 6451.953 (R 2 ), 6455.191 (R 3 ), 6459.463 (R 4 ), 
6465.369 (R 5 ) 

f remaining four heads are at positions 6774.2 (R 2 ), 6775.9 (R 3 ), 6779.6 (R4), 6 785.7 (R 5 ) dRosenll 19701) 

s wavelengths of the band heads extrapolated from spectroscopic constants of iHocking et all dl98d) : 6836.6 (R 2 ), 6839.9 (R 3 ), 6844.5 (R 4 ), 
6850.9 (R 5 ) 



gives good constraints on the extinction for regions with low- 
to-moderate reddening. For the position of V4332 Sgr, we get 
E B - V = 0.344 + 0.008. 



An independent constraint on the reddening can be obtained 
from the radio observations of the H 1 2 1 cm line. The integrated 
intensity of the HI line of 664.2 Kkm s can be converted to 
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21 



cm using a stan- 



a column density of A^(H I) = 1 .25 X 

dard conversion fac tor of 1.822 x 10 18 cm^ (K km s~')~' (e.g., 
Wils on et al.l 12009). The value of the redde ning can be found 
from a relation given by iBohlin et al.l (Il978l) . By neglecting a 
contribution from molecular hydrogen (our observations in the 
CO rotational transition, see Sect. 12.31 show no molecular cloud 
in the direction of V4332 Sgr) we get E B - V = 0.22 ± 0.04. In the 



present paper we adopt Eb-v - 0.32 as a reasonable compromise 
between the above estimates. 



The distance to V4332 Sgr is unknown . [Martini et al.1(ll999l) 
estimated that the object is at ~ 300 pc. Assumin g that the pro- 
genitor of V4332 Sgr was a main sequence star, iTvlenda et al.l 
( 2005) estimated a distance of ~ 1.8 kpc. A lower limit to the dis- 
tance can be obtained from the observed radial velocities of the 
interstellar absorption lines. A mean heliocentric radial velocity 
derived from the values listed in Table [3] is -5.5 km s~'. This 
corresponds to Vlsr — 6.5 km s~'. After adopting a Galactic ro- 
tation curve of Brand & Blitzj (1 19931) . we find that V4332 Sgr is 
at a distance k. 1.0 k pc. This excludes the low distance value of 
iMartiniet all (119991) . 
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As noted in iTvlenda et al.1 (l2005h . the radial velocity of 
V4332 Sgr shows that the object does not follow the Galactic 
rotation. The present analyses of the photospheric spectrum, 
atomic emission lines, and molecular bands give Vf, - -56 + 
1610ns- 1 (Sect, gjj, -65 ± Vkms" 1 (Table Q), and -75 ± 
10km s~' (Tabled, res pectively. These valu es are less negative 
than those estimated in Mart ini et ail fl999) (-180 km s _1 ) and 
Tvle nda et all (120051) (-160 km s _1 ), but are still noticeably dif- 
ferent from what is expected from the Galactic rotation in the 
direction of V4332 Sgr (k, -12 km s , heliocentric). Thus there 
is no way to estimate the kinematic distance of the object. 

In the present paper, we adopt, i f required, a distance of 
1.8 kpc (following lTvlenda et al.l l2005). This is a very uncertain 
value, but the principal conclusions drawn in the present paper 
are independent of the distance. 



5. Analysis of the molecular bands in emission 

The main question that arises when analysing the observed 
molecular bands in emission concerns mechanism(s) that could 
produce these features. One should consider collisional excita- 
tion and/or radiative pumping of excited states of the molecules 
as a primary source of the emissions. As is shown below, the 
observed band profiles can only be reproduced if the rotational 
temperature is assumed to be very low, typically ~120 K. This 
implies that the kinetic temperature must also be close to it in 
the ambient medium. In these conditions there is no way to sig- 
nificantly populate the observed upper electronic states by colli- 
sions, because they would require temperatures of a few thou- 
sand K, meaning that the only possibility that remains is the 
radiative pumping. Molecules at ground states absorb radiation 
photons capable of exciting them to higher electronic levels. 
Radiative decays from them produce the observed emissions. 
The analysis done in this section assumes this mechanism. 

For each molecule considered in the following subsections, 
we have constructed a molecular model including rotational lev- 
els of the ground and excited electronic states. Molecular data 
for the TiO bands were taken from Schwenke ( 1998). The source 
referen ces for VO and ScO are the same as in Kaminski et al. 
(2009). Generally, our models of the molecular structure in- 
cluded rotational levels up to the quantum number J — 65. 

The population of the rotational levels of the ground state, 
including spin substates, is assumed to follow the Boltzman law 
at a given rotational temperature. The upper electronic state (to- 
gether with its rotational levels) is assumed to be populated by 
radiative transitions from the ground state. The source of excit- 
ing radiation has usually been approximated by a blackbody of 
3200 K (expected effective temperature of an M5-6III star). In 
the case of the TiO bands, a synthetic spectrum from a model 
stellar atmosp here of T e ff=3200K a nd logg = 0.0 has been 
calculated (see Kaminski et al. 2009, for details). Atmospheric 
models of Hauschi ld et al] ([1999) were used for this purpose. 
Level populations under the above co nditions were solved with 
the code RADEX (Ivan der Tak et al.ll2007h . Assuming radiative 
decays of the excited electronic state, a synthetic emission can 
be calulated from the derived populations. Due to the noise of 
the observed spectrum, both the observed and calculated spec- 
tra have been smoothed out (typically with an FWHM of 2.4 A) 
before comparing them. 



5.1. TiO 

Numerous bands of titanium oxide are observed in our spec- 
trum. Usually narrow emission components at the positions of 
the molecular band-heads are superimposed on broad and deep 
absorption bands of the photospheric spectrum. The most promi- 
nent examples are three emission components of the TiO y (0,0) 
band at 7052, 7086, and 7129 A. We made a model analysis of 
these three emission components of the TiO y (0,0) band. The 
results are presented in Fig.|4] 

As can be seen from the figure, the emission components 
overlay a broad absorption photospheric component of the same 
band with a relatively sharp blue edge at the position of the F3- 
F3 emission component at ~7050 A. Therefore, for the purpose 
of the present analysis of the emission profile, we subtracted a 
synthetic stellar spectrum, calculated with T e ff - 3200 K and 
logg = 0.0, from the observed one (upper part of Fig. 0), The 
result is shown in the lower part of Fig. [4] 

The very first conclusion from the model analysis is that reli- 
able fits to the individual emission components can be obtained 
only if the rotational temperature is assumed to be 120±20K. 
For an increasing temperature, the red tail of each component 
becomes more pronounced and more extended toward longer 
wavelengths. 

Another conclusion can be drawn from the observed flux ra- 
tio, which is 2:1:3 for the F3-F3, F2-F2, and Fj-Fj components, 
respectively. Since the upper substates of each component are 
connected with the lower substates only by main branches (satel- 
lite ones are very weak), it can be assumed that all the three emis- 
sion components are uncoupled transitions. Then, differences in 
the fluxes can be caused by different populations of the lower 
substates and/or by differences in the flux pumping the upper 
levels. Assuming that the lower levels are populated according 
to the Boltzman distribution with a temperature of 120K, we 
find that the ratio of the Fi— Fj and F2-F2 fluxes is consistent 
with the optically thin limit. Effectively, the lowest rotational 
levels of the F3 substate (with energies from 400 to 520 crrT 1 ) 
are about 2 times less populated than the lowest rotational levels 
of the F2 substate, and 4 times less than those of the lowest Fj 
substate. If the pumping flux is similar for the three transitions, 
the expected value of the above flux ratio is therefore 1 :2:4. The 
observed F3-F3 component is thus significantly enhanced. An 
assumption of optically thick emissions would only flatter this 
ratio. The enhanced emission of the F3-F3 component could be 
explained if the flux illuminating the emitting area in this tran- 
sition were significantly higher than the fluxes pumping the re- 
maining two. As can be seen from Fig. [4] (upper part) the F3-F3 
component (at ~ 7050 A) is situated on the strong blue edge of 
the photospheric absorption band. Thus, if the emitting region 
escapes from the source of the photospheric radiation, the edge 
of the photospheric band is redshifted and a significant enhance- 
ment of the pumping radiation is then available for the F3-F3 
component compared to the F1-F1 and F2-F2 components. Our 
modelling shows that the observed flux ratio can be satisfactorily 
reproduced if the photospheric spectrum, as seen by the emitting 
region, is redshifted by at least 25kms~'. The resulting emis- 
sion profiles, calculated for a relative velocity between the pho- 
tospheric source and the emitting gas of 40 km s" 1 , are shown in 
the lower part of Fig. [4] 

5.2. ScO and YO 

In the observed spectrum, both ScO A 2 n-X 2 £ + (0,0) and YO 
A 2 II-X 2 Z + (0,0) bands are present in emission. The rotational 
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VO B 4 n-X 4 E (0,0) 





y V VY 



7100 
Wavelength (A) 



7900 7950 
Wavelength (A) 



Fig. 4. Top panel: the observed spectrum in the vicinity of the 
TiO y (0,0) band (black line) and a synthetic model spectrum for 
r e ff = 3200 K and \ogg = 0.0 (red line). Bottom panel: result 
of subtraction of the synthetic stellar flux from the observed one 
(black line) and the model emission calculated as explained in 
the text (blue line). 



Fig. 6. The observed spectrum near the emission band of VO 
B 3 II-X 4 5L (0,0) (in black). A spectrum of M7 III (HD 18242) 
star is overplotted in blue. The sum of the M7 III spectrum and 
of the synthesized rotational contour is shown in red. See text for 
details of the model. 




6060 
Wavelength (A) 



Fig. 5. The observed spectrum around the ScO A 2 Z-X 2 E + band 
(black line). The probable identification of the CrO B 5 II-X 5 n 
(0,0) and (1,1) bands is also marked. A model emission profile 
of ScO is shown in red. 



contours of ScO sub-bands, A 2 n 3/2 -X 2 £ + and A 2 n I/2 -X 2 £ + , 
are shown in Fig. [5] A grid of theoretical profiles have been syn- 
thesized for a range of rotational temperatures. The best fit to 
the observed spectra is obtained for a rotational temperature of 
120+20K and a radial velocity of -TSilOkms" 1 (Fig.0. 

Emission of yttrium oxide A 2 ni /2 -X 2 £ + (0,0) at 6133 A is 
rather weak and barely recognizable. The second component 
A 2 n 3/ 2-X 2 2 + (0,0) of YO is not seen in the spectrum. 

5.3. VO 

Vanadium oxide B 4 n-X 4 5T (in short B-X) (0,0) and (1,0) bands 
are seen as prominent emission features (see Fig.Q]and Table|2]i. 
The most intense B-X (0,0) band is shown in detail in Figure 



[6] The band consists of four sub-bands formed by the main 
branches. A model fit to the data suggests a rotational temper- 
ature of ~100 K, in agreement with the rotational temperature 
derived above for the TiO bands. 

The observed spectrum around the VO B-X (0,0) band has a 
relatively high S/N. Comparing this part of the spectrum with 
the UVES spectrum of HD 18242 (M7 III), we can estimate 
the radial velocity of the photospheric spectrum. From the sharp 
absorption feature (probably due to VO in the photosphere) at 
7960.4 A, as well as the nearby TiO band head at 7818.6 A, we 
estimate the heliocentric stellar velocity to be -60±20 kms~'. 

A final fit was obtained by summing the predicted rotational 
countour of the VO B-X band and the spectrum of HD 18242. 
The best fit (shown in red in Fig. [6} was obtained for a velocity 
of emitting gas of -75+5kms~'. The observed spectrum also 
shows the VO C 4 5T-X 4 5r (1,0) and (0,0) bands around 5474 A 
and 5736 A, respectively. They appear as narrow emission fea- 
tures, which are too weak and too noisy to be analysed in more 
detail. 

5.4. CrO 

After having identified the atomic and molecular emission fea- 
tures in the observed spectrum, we are left with a few broad 
(~ 20 A) emission-like features in the spectral range between 
5800 and 6800 A. We postulate that they come from the elec- 
tronic band system B 5 n-X 5 n of CrO observed in emission, as 
shown in Fig Q] Details of the proposed identification can be 
found in Table [2] 

We have identified five vibronic bands of CrO, i.e. (1,0), 
(0,0), (0,1), (0,2), and (2,3), shown in Figs. Q and |5] The ex- 
pected positions of the (1,2) and (1,1) bands are also marked in 
the figures. The broad emission profiles are due to overlapping 
of five sub-bands formed by main branches. 

Relative intensities of the emission bands with a common 
upper vibrational state can be compared with theoretical pre- 
dictions. For this purpose we used Franck-Condon factors com- 
puted bv iReddv et al] (Il998l) . We find that the predicted ratios 
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Wavelength (A) 
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O-0 B 5 n-x 5 n (0,1) 
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Fig. 7. Detailed view of CrO B 5 E[-X 5 n bands in emission. The 
observed spectra are plotted with a black line, while theoretical 
contours of the CrO bands are shown as a blue line. An M7 III 
spectrum is displayed in magenta in the middle and bottom pan- 
els. A synthetic spectrum is shown in green in the top panel. 



are fairly consistent with the observed fluxes. The ratios of the 
emission fluxes of the (1,0):(1,1):(1:2) bands are predicted to 
be 3.0:0.1:2.0, while the observed ones are 7:0:^1, respectively. 
The observed ratios of the (0,0):(0,1):(0:2) bands are 3.5:4.8:3.4 
compared with the theoretical predictions of 3:4:2, respectively. 
We synthesized the expected contours of the CrO emission 
feature s. A line list from the laboratory work of Hocking et al. 
(1980) was used. Rotational intensities were computed con- 
sistently f rom diagonalization of Hamiltonians defined by 
Hocking et al. ( 1980) and using the standard approach (see e.g. 
Hougen 2001). The computed strenghts were checked to re- 
produc e analytical formulas of rotational intensities of Kovacs 
(1969). The molecular structure included in the model is rather 
complex and consists of the four lowest vibrational levels of the 
ground state and the three lowest vibrational levels of the upper 
electronic state. Each vibrational level is composed of five sub- 
states and each substate is composed of lower rotational levels 
up to J - 65, in order to include levels with energies up to 2400 
cirT 1 . All the radiative transitions between the listed levels are 
included. 



To reproduce the observed shapes of the rotational contours, 
a rotational temperature of the ground state of 450 K was as- 
sumed. This relatively high rotational temperature (compared 
to that derived from the above analysis of other molecules) is 
needed to explain formation of all the five band heads seen in 
the observed spectrum, as well as their trapezoidal shapes. The 
exciting radiation was simulated by a blackbody of 3200 K. The 



obtained theoretical contours are shown in Fig. [7] They have 
been shifted by -70 km s to fit the observed features. 

Because of the reasonable consistency between the observed 
and predicted wavelength positions, rotational contours (as ex- 
pected from environment of low kinetic temperature), intensity 
ratios, and the presence of all the expected bands, we consider 
the identification of the CrO B-X system in the spectrum of 
V4332 Sgr as reliable. CrO was sai d to be seen in the atmosphere 
of p Pegasi (spectral class M2) by Davies ( 1947) on the basis of 
band-head mesurements. Later, on this finding was questioned 
by Jorgensen (1996). Our identification is the first detection of 
CrO in emission in any astrophysical environment. 



6. Origin of the emission features 

Two important conclusions have been drawn from our model 
analysis of the observed molecular bands in emission in Sect. [5] 
First, these features are produced by radiative pumping. Second, 
the emitting gas is escaping from the source of the pumping ra- 
diation (Sect. ED. 

If the molecular emission bands are formed by radiative 
pumping then it is natural to suppose that the atomic emission 
lines are also produced in the same mechanism. All the atomic 
lines observed in our spectrum in emission are from resonance 
transitions. In this case it is difficult, if not impossible, to avoid 
radiative pumping in the atomic lines, if this mechanism is effec- 
tive for the molecular bands; both photoabsorption cross sections 
and the expected abundances are larger for elements like Na and 
K than for molecules like TiO, VO, or ScO. 

As can be seen from TableQ] the observed ratio of the Na I D 
lines (5890/5896 A) is close to 1.0. The same ratio for the KI 
lines (7665/7699 A) is ~1.2. In an optically thin case both ratios 
are expected to be 2.0. The fact that the line ratio in the Nal and 
K I doublets is observed to be close to 1 .0, can be very easily and 
naturally explained within the hypothesis of radiative pumping. 
The line intensity in this mechanism is limited by the flux of the 
pumping spectrum. Since the number of photons available for 
each line of a given doublet is practically the same in the pump- 
ing spectrum (assumed to be stellar-like), already for a moderate 
thickness (of a few) in the lines, the line ratio must approach 1 .0. 

In terms of intensities, the K I doublet dominates the spec- 
trum in the observed spectral range. It is ~25 times stronger 
that the Na I doublet. The high K I/Na I flux ratio is rather un- 
usual. In most circumstellar environments, when both doublets 
are observed, the Nal lines are stronger thanks to the ~16 times 
higher abundance. (The transition probabilities for both doublets 
are very close.) The relative KI and Nal line strengths observed 
in V4332 Sgr can be easily explained, if they are assumed to be 
pumped by radiation. As discussed above, the observed line ra- 
tios indicate that the lines absorb almost all the available pump- 
ing radiation. Therefore the observed K I/Na I flux ratio is ex- 
pected to reflect the spectral energy distribution of the underly- 
ing photospheric source. The flux ratio of the observed photo- 
spheric continuum near the KI and Nal doublets can be esti- 
mated from the photometry in the Iq and V bands, respectively 
(effective wavelengths of these photometric bands are close to 
the wavelengths of the doublets). By adopting Ic = 15.1 and 
V = 19.7 (see Sect. 13 and using the standard Vega flux cali- 
bration, one gets ^(7680 A)/AF A (5890A) =* 28. Another way 
to estimate the above ratio is to measure the continuum from 
our Subaru spectrum. Near the Na I doublet, the observed con- 
tiunnum is very faint but it can be estimated as F / i(5890A) =* 



1 X 10 "ergs" 



In the vicinity of the K I doublet, we 
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measure F^(7680A) - 1.7 x 10~ 16 erg s -1 cm -2 A" 1 . This results 
in AF A (l 680 A) /AF A (5890 A) =* 22. Both estimates are in excel- 
lent agreement with the observed flux ratio of the K I and Na I 
doublets, confirming our hypothesis that the emission lines ob- 
served in the spectrum of V4332 Sgr are indeed pumped by the 
photospheric radiation of the central source. 

As can be seen from Table Q] the EWs of the Ca I and Rb I 
lines are almost two orders of magnitude lower that those of the 
K I lines. This implies that the former lines are optically thin; i.e., 
they absorb only a small part of the available photospheric pho- 
tons. This statement is confirmed by the expected low abundance 
of Rb and the low transition probability in the Ca I (intercom- 
bination) line. The difference in the optical thickness is proba- 
bly also the reason for the observed differences in the profiles, 
namely that the Cal and Rbl lines are significantly narrower 
compared to those of KI and Nal (see column 6 in Table[TJ. 

We have shown that the spectral shape of the observed pho- 
tospheric spectrum of M5-6 III can explain both the shape of the 
TiO y (0,0) band (Sect. ED and the flux ratio of the KI/Nal 
doublets. However, in the observed spectrum we do not see any 
absorption features, which could be interpreted as from photons 
absorbed in the pumping mechanism. This shows that we do not 
see any absorbing matter along the line of sight and implies that 
the matter seen in the emission features is not distributed spher- 
ically around the continuum source. The observed continuum 
is also much too faint to be able to pump emission features as 
strong as observed. As can be seen from Fig.Q] peak intensities 
in some molecular features and particulary in the atomic lines 
are by large factors higher than the intensity of the underlying 
continuum. The observed EWs of the K I lines (see Table [TJ 
imply that if they had been excited by the observed continuum 
in their vicinity, they would have had to absorb all the contin- 
uum photons from a spectral region as large as ~550 A. This 
would be possible only if the absorbing matter had been moving 
with a velocity gradient, between slowest and fastest regions, 
as large as ~20000 km s , which is two orders of magnitude 
highe r than the velocitie s observed in V4332 Sgr during its erup- 
tion ((Martini et al. 1999). Therefore we conclude that the photo- 
spheric continuum responsible for pumping the observed emis- 
sion features is much stronger than the observed one. The above 
discussion of the KI lines suggests that we observe ^1% of the 
flux emitted by the central object. In Sect. [7] we further argue 
that the central object in V4332 Sgr is mostly hidden to us by an 
opaque disc seen almost edge-on. 



7. Spectral energy distribution 

The SEP o f V4 332 Sgr was already investigated by 
iTvlenda et al.l d2005f) . It was based on photometric data ob- 
tained in 2003 and covering photometric bands from B to M. 
Tvlendaetal. (2005) point out that apart from the stellar-like 
component dominating in the optical and near-IR, the object dis- 
played a large IR-excess seen in the KLM bands. When inter- 
preted as a single blackbody the IR component had an effective 
temperature of ~750 K and a luminosity ~15 times higher than 
that of the optical co mponent. Another possibility mentioned in 
ITvlenda et alJ d2005h was an accretion disc ~25 times more lu- 
minous than the central star. 

Figure[8]displays an SED of V4332 Sgr derived from various 
observations obtained between April 2005 and May 2009. The 
source of the data are as follows: 



- Optical photometric measurements were taken from the web 
page of V. GoranskiJ_|. These are combined data from differ- 
ent observatories. We interpolated the BVRcIc magnitutes 
for the date of the Subaru/HDS observations and adopted 
B=21.4, V=19.7, R c =n.9, / c =15.1 mag. The error bars 
shown in Fig. [8]for the optical magnitudes represent a scat- 
ter in the photometric measurements around the date of the 
Subaru/HDS observations. 

- The JHKs data are from the observations reported on in 
Sect.O 

- IR measurements in the LMN„ bands were taken from 
iLvnch et all (12006). and they correspond to observations car- 
ried out on 7 August 2006. 

- Far-IR data were taken from iBaneriee et al.l d2007l) . These 
are Spitzer/MYPS fluxes at 24 //m, 70 pm, and 160 pm aver- 
aged from measurements obtained on 15 October 2005 and 
2 November 2006. 

- AKARI measurmements were taken form the first releases 
of the FIS Bright Source Catalogue ( Yamamura et al.ll2010h 
and IRC Point Source Catalogue ( Kataza et al.ll2010l) . which 
are based on observations from the AKARI All-Sky Survey 
performed in 2006-2007. For the AKARI Infrared Camera 
(IRC), the source was detected only in the 9 pm band 
with flux density of 507.7 mJy. Reliable measurements with 
the Far-Infrared Surveyor instrument (FIS) are available 
only for the 65 pm and 90 pm bands with density fluxes 
of ^65 = 1.79 Jy and 7^90=1.13 Jy. There are also data for 
the FIS bands at 140 and 160 pm (Fi 40 =86.9 mJy and 
Fi6o=29.9 mJy), but the measurements are flagged as un- 
reliable (too weak a source at those bands). The total (ab- 
solute and relative) calibration uncertainty of the measure- 
ments in all the AKARI bands is 20% (3cr; I Yamamura et all 
l20ToHKataza et alJl20Toh . 

- In Fig. [8] we also display spectra of V4332 Sgr in wide spec- 
tral ranges. The optical spectra are the Subaru/HDS spec- 
trum reported i n this paper and a synth etic spectrum from the 
MARCS grid dGustafsson et alJl2008l) with T eff = 3200 K, 
solar metalicity, and log g - 0.0; the last was rescaled to fit 
the fluxes of the Subaru/HDS data. The Spitzer spectra are 
also shown. They were acquired on 18 April 2005 and 27 
September 2005. 

Vega fluxes were used as zero points when converting 
magnitudes to flux-density units. The optical spectra and the 
BVRqIcJHKs data points shown in Fig. [8] are corrected for in- 
terstellar reddening with Eg-v - 0.32 and Ry = 3.1. 

In the following discussion we assume that the object did not 
vary significantly between 2005 and 2009. This is not entirely 
justfied since, as noted in Sect.Q] the object faded in the optical 
in 2006; however, it seems that the overall SED has remained 
more or less unchanged over a few recent years. When compared 
with the photometric data in 2003 (Tvlenda et al. 2005), the ob- 
ject is now ~1.5 mag fainte r in BVR but ~1.0 mag brighter in 
JHK. iBanerjee et al.l d2007l) note that the object remained con- 
stant over a year in their far-IR observations in 2005-2006. 

Similar to lTvlenda etal] d2005l) . the SED in Fig.[8]displays a 
huge IR-excess. However, thanks to the far-IR data, we can now 
conclude that it cannot be represented by a single temperature 
component. We have attempted to reproduce it in Fig. [8] with 
a sum of two blackbodies of 950 K and 200 K, but it would 
be more reasonable to conclude that this is a multi-temperature 
component with temperatures ranging between these two values. 
This result, combined with the fact that the IR component has 



http://jet.sao.ru/~goray/v4332sgr.ne3 
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Fig. 8. Spectral energy distribution for V4332 Sgr. Observational 
points are taken from different sources as listed in the text. Our 
attempt to fit the observed SED is shown with a dashed line. 



remained stable since at least 2003, strongly suggests that we 
see a multi-temperature Keplerian disc of matter surrounding the 
central object resembling an M5-6 giant. 

In the simple fit to the observed points in Fig. |8]the two IR 
components, i.e. blackbodies of 950 K and 200 K, are ~30 and 
~20 times more luminous, respectively, than the stellar compo- 
nent of T e ff = 3200 K. Assuming spherical geometry of the com- 
ponents, the two infrared components would have radii ~65 and 
~ 1 200 times larger than the stellar one. If interpreted with a disc- 
like spectrum, the disc would have an apparent luminosity ~50 
times higher than the central star. 

In principle, one can consider two possible sources of energy 
radiated away by a disc, i.e., reproduction of the central star ra- 
diation absorbed by the disc matter (so-called passive disc) or 
gravitational energy dissipation by viscous forces transferring 
matter from outer disc layers to inner ones and eventually to 
the star surface (accretion disc). In both cases, the central star 
is expected to be at least as bright as the disc. In the first case 
(passive disc), for obvious geometrical reasons the disc can ab- 
sorb, hence reradiate, only a part of the stellar luminosity. In the 
second case (accretion disc), the accreted matter is expected to 
release a comparable energy to what is dissipated in the disc, at 
the star surface (so-called disc boundary layer), where the matter 
has to slow down from the Keplerian velocity to the velocity of 
stellar rotation (usually much lower than Keplerian). In the case 
of a giant or a supergiant, the boundary-layer braking is likely 
to take place below the photosphere. In this case, accretion is 
expected to add a luminosity to the intrinsic stellar luminosity, 
similar to that of the disc. 

The only way to reconcile that the central star in V4332 Sgr 
is observed to be -50 times less luminous that the disc is to 
assume that the disc is observed almost edge-on. In this case 
the disc would block most of the stellar radiation in the direc- 
tion to the observer. The geometrical aspect (inclination close to 
90°) means that the observed lumnosity of the disc is also ex- 
pected to be significantly lower than the true one. Therefore we 
can conclude that the stellar component observed in the optical 
only accounts for %\% of the luminosity of the central object in 
V4332 Sgr. This conclusion nicely fits the one drawn from our 
analysis of the emission features in Sect. [6] 



8. Discussion 

Our analysis done in Sects. [6] and Q has led us to conclude that 
V4332 Sgr is now observed as consisting of a central object, re- 
sembling an M5-6 giant, surrounded by a disc of circumstellar 
matter. The inclination angle of the disc is close to 90°, so that 
the disc matter obscures the central star, and we can only observe 
a tiny part of its radiation. The situation thus resembles what is 
observed in IRAS 18059-3211, known as Gomez's Hamburger 
dRuiz et al.ll 19871 IWood et al.ll2008l) . whose global SED is simi- 
lar to V4332 Sgr. The central star in Gomez's Hamburger is not 
directly seen, and the observed optical spectrum is entirely due 
to scattering of the central star radiation on dust grains in the 
outer edge of the disc. It is very probable that the same situation 
is also observed in V4332 Sgr. Polarimetric observations in the 
optical would allow verification of this hypothesis. 

Molecular bands in emission have been observed in the op- 
tical spectrum of the luminous red supergiant VY CMa (Herbig 
1 1 974t IWal lerstein & Gonzalez 2001 , and references therein), as 
well as of U Equ (Barnbaum et al. 



1996). In both cases it has 



been suggested that the objects are surrounded by massive cir- 
cumstellar discs, which obscure the central stars along the line of 
sight (Herbig 1970; Kastner & Weintraub 1998; Barnbaum et al. 
1996), i.e. similar to what we propose in the present paper for 
V4332 Sgr. The existence of a disc obscuring the central star 
thus seems to be a condition favouring observations of molecu- 
lar bands in emission. Indeed, if we increased the photospheric 
spectrum in V4332 Sgr by a factor of 100 (as suggested in 
Sects. |6]and|7] if the central star were unobscured) remaining the 
fluxes of the emission features unchanged, then we would have 
practically no chance to detect any emission molecular band. 
The only visible emission features would then be those of KI 
andNal. 

As shown in Sect. 13.11 the spectral type of the observed pho- 
tospheric spectrum agrees with the BVR C I C photometry corrected 
for the interstellar extinction. Thus the scattering responsible for 
the observed photospheric spectrum has no significant effect on 
its global shape. This can be understood if the observed spectrum 
is mainly produced by forward scattering on the outer disc rim 
lying in front of the central star, as we suggest for V4332 Sgr. In 
this case scattering is dominated by large grains (a ^ 0.3 jjm), 
which have a high value of g =< cos 6 >, i.e. scatter mainly in 
forward directions. These grains also have Q SC!i th at is practically 
independent of wavelength in the optical (e.g. lLoar & Draind 
I1993 )F1 

As discussed in Sect. [6] the emission features observed in 
V4332 Sgr are most likely pumped by radiation from the cen- 
tral star. One may ask where the matter seen in the emission 
features can be situated. One possibility is that this could be 
matter in the atmosphere of the disc. However, the profiles of 
the emission lines observed in V4332 Sgr are single-peaked (see 
Fig. [3}, while it is widely believed that emission features ob- 
served from an edge-on disc should be double-peaked (due to 
the disc rotation). Although this is true in most cases, there are 
situations where disc e mission lines can be observed as single- 
peaked. Simulations of [Murray & Chiang! (1 19961 1 19971) show 
that this happens when the disc is seen at a high inclination an- 
gle and possesses a wind, which is optically thick in the lines. 
This mechanism could explain the profiles of the Na I and K I 
lines, since, as discussed in Sect. [6] these lines are most likely 
optically thick. However, as discussed in the same section, the 
Ca I and Rb I lines are optically thin, and yet they are single- 

8 see also http://www.astro.princeton.edu/~draine/dust/dust.html 
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peaked (see Fig. [3). Another problem arises from the observed 
line wings extending up to ~ ±250 km s from the line cen- 
tre (see Fig. [3}. If formed in a high-inclination disc, they would 
measure the rotation velocity of the innermost observable re- 
gions of th e disc. Adopting a d istance of V4332 Sgr of ~ 1.8 kpc 
(following lTvlenda et al.ll20Q5l who obtained this value after as- 
suming that the progenitor was a solar type star), our fit to the op- 
tical photometry (mentioned in Sect. 13.11 ) results in an effective 
radiu s of ~5 R (similarly to the one derived in Tvlend aet al.l 
2005, for 2003). However, as we now know, the true central star 
is probably two orders of magnitude brighter than observed and 
thus has a radius that is an order of magnitude larger than the 
above estimate. For a solar-mass star and a radius of ~50 R one 
obtains a Keplerian velocity of ~60 km s _1 at the stellar surface. 
Since the central star is obscured and thus the disc regions up 
to a few stellar radii are also unseen, there seems to be no way 
to explain the line wings extending up to ~ ±250 km s within 
the idea of the emission lines arising from the disc. A higher 
stellar mass would make the situation even worse. For instance, 
a 3 M progenitor would be almost two orders of magnitude 
brighter and about 10 times more distant, and the outburst rem- 
nant would be 10 times more extended, lowering the Keplerian 
velocity at its surface by a factor of 2. 

It thus seems that we are only left with a scenario, in which 
the emission features are formed in matter lost by the central ob- 
ject. This idea is in accord with our findings in Sect. l5.1l that the 
observed ratio of the emission components of the TiO y (0,0) 
band implies an outflow velocity of the emitting matter of at 
least 25 km s . The origin of the outflowing matter could be 
matter ejected in the past, e.g. during the 1994 outburst, or a 
wind flowing out at present from the object, or both. The low 
rotational temperature found from the molecular bands in emis- 
sion (Sect. |5J suggests that these features are formed in a cold 
medium, thus at significant distances from the central object. 
This seems to favour the former possibility of matter ejected in 
the past. The observed widths of the atomic lines at zero inten- 
sity (see Fig.0 imply an outflow velocity of up to ~ 250 km s _1 . 
The line profiles suggest that the outflow is concentrated along 
an axis nearly perpendicular to the line of sight. This could be 
the axis of the disc. Detailed modelling of the line profiles would 
be required to better constrain the structure and nature of the out- 
flow. 
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